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Abstract 

Quantum walks, both discrete (coined) and continuous time, form the basis of several quan- 
tum algorithms and have been used to model processes such as transport in spin chains and 
quantum chemistry. The enhanced spreading and mixing properties of quantum walks compared 
with their classical counterparts have been well-studied on regular structures and also shown 
to be sensitive to defects and imperfections in the lattice. As a simple example of a disordered 
system, we consider percolation lattices, in which edges or sites are randomly missing, interrupt- 
ing the progress of the quantum walk. We use numerical simulation to study the properties of 
coined quantum walks on these percolation lattices in one and two dimensions. In one dimension 
(the line) we introduce a simple notion of quantum tunneling and determine how this affects the 
properties of the quantum walk as it spreads. On two-dimensional percolation lattices, we show 
how the spreading rate varies from linear in the number of steps down to zero, as the percolation 
probability decreases towards the critical point. This provides an example of fractional scaling 
in quantum walk dynamics. 
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1 Introduction 



Quantum versions of random walks have been extensively studied since they were introduced in the 
context of quantum algorithms A classical random walk is essentially a diffusion process, in 

which spreading from the initial configuration occurs as a series of random steps. In a quantum walk, 
quantum coherences replace the diffusion, which can lead to both faster and slower spreading from 
the starting point of the quantum walk. Faster spreading provides the speed up in man y quantu m 
walk algorithms. Exponential speed up has been proved for certain transport problems: Kempe 3; 
[H] showed h ow a quanturn walk can cross a hypercube much more efficiently than a classical random 
walk, while Childs et al. 0] have produced a scheme for a continuous time quantum walk that 
can find its way across a particular "glued trees" graph exponentially faster than any classical 
algorithm. On the other hand, the phenomenon of Anderson localizati on of quantura particles is 
we ll-known, and ha s been related to continuous-time quantum walks by iKeating et al.l 7|] . Studies 
by Krovi and Brun 0, [13] extend this to coined quantum wa lks, and highlig ht the importance 



of symmetry in determining the properties of quantum walks. IStefaiiak et al.l 111 . Il2l | provide a 



way to characterize the localization of unbiased quantum walks on regular lattices, showing how it 
depends on the topology of the lattice, the coin operator and chosen initial state. 

Localization can also be exploited, as in the quantum walk search algorithm, wh ere a marked 



state perturbs the uniform distribution causing the quantum walk to converge onto it. IShenvi et al, 



[lit proved that basic quantum walk searching is quadratically faster, the maximum speed up that 
can be achieved for unordered searching IJ]. This quadratic improvement is only achievable for a 
well-connected database though, on a two-dimensional grid it is still open whether the full quadratic 
spe ed up is po ssible. For a recent review of quantum walk searching and related open problems. 



see 



Santhal (ij. 



Quantum walks also provide useful models of physical phenomena such as quantum state trans- 
fer in spin chains 16|], or energy transport in biomolecules [17[. Here the goal is to transfer a 
quantum state or an excitation, with high fidelity being more important than high speed. Cer- 
tain combinations of parameters can achieve perfect state transfer (reviewed in [18|), while a small 



amount of noise can produce robust near-perfect transfer in more disordered systems 13] . There 
is no general formula for predicting when a quantum walk will be more efficient than its classical 
counterpart, and the wide-ranging studies to date (e.g., [l^) only scratch the surface of the possible 
configurations and parameters available. 

In this work, we are motivated by using quantum walks as models of transport phenomena, so 
we are interested in how fast they spread from a localized starting state. In order to have a simple, 
but flexible, concrete setting for our work, we concentrate on the behaviour of quantum walks 
on percolation lattices, i.e., Cartesian lattices with vertices or edges randomly missing, that thus 
limit the paths the quantum walker can take. Edge (bond) and site (vertex) percolation lattices 
are defined by the probability p of an edge or site respectively being present in the structure. 
In two or more dimensions, as the proportion of missing edges or sites decreases, around some 
critical probability pc the lattice changes from a collection of disconnected sections to almost all 
connected. Moreover, there will in general be a path from one side of the lattice to the opposite 
side for p > pc, an obvious qualitative change in the transport properties of the lattice. The two 
types of lattice (with either site or bond percolation) have similar properties, but different critical 
parameters. They are used to model the behaviour of phenomena as diverse as the spread of 
contagious diseases, the propagation of forest fires, and the movement of oil deposits in porous 
rocks. They have been studied in great detail and, where analytical solutions are not available, 
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efficient numerical methods have been develop ed 20. |21|. For an ac cessible introduction to the 
properties and uses of percolation lattices, see StaufFer and Aharony [23 |. For this work, we are 
simply using the percolation lattices as a substrate for the quantum walks. One of the most basic 
questions we cai i ask is how the quantum walk behaves as the density of edges or sites is varied. 
In related work, IXu and Liul 231 ] study quantum walks on Erdos-Renyi networks, which also have 
randomly placed edges similar to percolation lattices, but without the underlying regular structure 
of a lattice. The limiting case is thus a walk on the complete graph, rather than the more familiar 
walk on a line or grid. 

Like classical random walks, quantum walks come in both discrete time ^MMM, and con- 
tinuous time {l| versions. Here we investigate only discrete time (coined) quantum walk dynamics, 
though we expect the continuous time dynamics to be qualitatively similar in many respects. We 
first review the properties of a discrete quantum walk on the line. For the first half of the paper, 
we work mainly with edge (bond) percolation lattices. In one spatial dimension (the line), just a 
single missing edge is enough to halt the progress of both quantum walks and random walks. So we 
explore two less restrictive variations on this theme: dynamic gaps, where the missing edges change 
with each time step |26l ]: and tunneling, a simple model of the quantum phenomenon whereby a 
quantum particle is able to pass through a narrow barrier with some small probability. There is 
also no (non-trivial) phase transition in one dimension for the percolation lattice, effectively = 1- 

We follow these one-dimensional examples with studies on two dimensional lattices. In general 
there is now more than one possible path between two points, allowing for non-trivial behaviour by 
the quantum walker, even with gaps fixed throughout the walk. Bond (edge) percolation in a two- 
dimensional Cartesian lattice has a phase transition at pc = 0.5. We also compare site percolation, 
where pc = 0.5927. . . for two dimensions 2^, 21]. We then discuss our results in the final section. 



2 Walk on the line 

A discrete quantum walk on a line can be defined in direct analogy with a classical random walk, 
where the walker repeatedly steps one unit forward or back based on the outcome of a random 
coin toss. The quantum walker has a quantum coin that will in general be in a superposition of 
"forward" and "back", and the quantum walker can step into a superposition of positions on the 
line, based on the state of the coin. However, the coin toss is not random, since pure quantum 
dynamics need to be unitary, but if measured, the outcome of the measurement would be random, 
like the classical coin. 

To fix our notation, we define a two-dimensional quantum coin with basis states ]c) with c E 
{±1}, and label the position basis states \x) by the integer location on the line x. A general state 
of the quantum walker can thus be written 

\'ip{x,t)) =^ac,x{t)\c,x), (1) 

c,x 

where ac,x{t) is the complex amplitude associated with the walker being at position x with the coin 
in state c at time step t, satisfying the normalization condition Ylcx l^c,x(i)l^ = 1- Here we have 
written \c,x) = \c) \x) for combined basis states. The evolution of the walk is governed by a 
coin operator C2 that acts on the quantum coin at each step of the walk. The simplest and most 
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commonly used coin operator is the Hadamard operator, 



,(Had) 



1 

7^ 



1 1 
1 -1 



(2) 



After "tossing" the coin with the coin operator, the particle moves to adjacent positions according 
the the coin state; this is expressed mathematically as a conditional shift operator 



S|-l,x) 
S|+l,x) 



1+1, X + 1) 



One step of the quantum walk is produced by the unitary operator U 
walk of t steps starting from an initial state IV'o) is thus 

mx,t)) = u'm- 



S(C (g) 



(3) 

A quantum 
(4) 



The position pr obability distribution of a quantum walk on a line, originally solved analytically 
by Ambainis et alj 0;[23]) is by now well-known, the time evolution over 100 steps with a Hadamard 
coin operator and initial state of -^(1+1,0) -|- i| — 1,0)) is shown in figure [H The double-peaked 
spreading is quite different to the binomial distribution of a classical random walk, and expands at 
a linear rate, giving a quadratic speed up over the ^/t spread of the binomial distribution. 




' on line thw 
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100 



Figure 1: Probability distribution for a quantum walk on the line over 100 steps showing the 
time evolution, using a Hadamard coin, equation ([2]), and a symmetric initial state (|-|-1,0) -|- 
i| — 1, 0))/^/2. Only positions with non-zero probability of occupation are shown, since odd positions 
are unoccupied at even time steps and vice versa. 
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Since the walk on the hne is amenable to analytic treatment, it has been exhaustively stud- 
ied by many authors o yer over the pas t ten years. Both path-counting and Fourier transform 



method s were given by Ambainis et al. 0], and path counting (path integrals) were further re- 



fined bylCarteret et alJ |28l |. A method using the algebra of the matrix operators was presented by 



Konno et al 



29l. 30ll a nd the tools of classical optics were adapted to quantum walks by Knight et al, 



3ll |. iRomanelli et alJ 13'Z\ separate the dynamics into Markovian and interference terms, which they 
relate to the dynamics of a kicked rotor. The coin operator may be generalised to any SU(2) oper- 
ator, 



T] 



(5) 



where < 6,(j) < ir are arbitrary phase angles and the coin bias is < r/ < 1. If we set rj = 0.5 and 
9 = (p = 0, we obtain the Hadamard coin operator of equation ([2]). The initial state of the coin can 
be similarly varied, 

1^(0)) = 7^1+1, 0) + e^^Vl-^'|-l>0), (6) 

where < 0' < tt is an arbitrary phase angle and the coin initial state bias < r/' < 1. As was 
shown bv lBach et al.l 33|, it is sufficient to specify one phase to obtain the full range of behaviour 
from the subsequent walk on the line starting from the origin. It is convenient for our purposes 
here to place the phase factor in the coin initial state {9'), and deal only with coin operators of the 
form 



(7) 



Varying the phase alone {rj = rj' = 0.5) produces biased walk distributions with a greater probability 
of the walker being found in the positive (negative) directions for 9' = 0(7r). Setting 9' = Tr/2 gives 
the unbiased walk shown in figure [TJ Varying the bias r/ in the initial state also controls how skew 
or symmetric the resulting walk turns out to be 3J, |35|. Varying the bias in the coin operator 
controls how fast the walk spreads. The extreme cases are t] = 1, for which the walker 
hops along the line without any reversing, in one or both directions as indicated by the initial 
coin state, and rj = 0, for which the walk oscillates back and forth between the initial state and 
nearest neighbouring positions, making no further progress along the line at all. There have also 
been exte nsive stu dies of decoherence applied to quantum walks, both analytic and numerical, see 
review by Kendon 3^, and references therein. 



Numerical simulation of a quantum walk is a straightforward evaluation of equation dH) for 
chosen parameters. In high level computational environments such as Matlab, it takes only a few 
lines of code. However, we wrote our own routines in C and C-|--|- to ensure they were efficient 
enough to to run up to at least 10,000 time steps for the walk on the line, and over 100 steps for 
the walk on a 2D lattice (the memory required scales linearly in the number of lattice sites). It was 
also necessary to repeat the simulations over a sufficient number of random percolation lattices to 
ensure reasonable statistics were obtained. The largest of our simulations took over a week on a 
SGHz processor to complete this averaging for 5,000 random lattices. This was in fact largely due 
to the time required to generate t he random numb ers for specifying the percolation lattice (one 
per walk) . We used routines from iPark and Milleii 37l | freely available online for generating the 
random numbers, although the quality of the random numbers is not particularly important for 
this type of simulation. 
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2.1 Dynamic gaps 



If edges are missing from the line, neither the classical nor the quantum walkers have any way to 
move beyond the gaps, and are thus constrained to the section of the line on which they start. 
A more interesting problem is the effect of gaps that change location after each time step. The 
gaps are then only a temporary barrier to the walkers' progress. In order to apply a quantum walk 
to a dynamically changing structure, the unitary evolution becomes dependent on both position 
and time. The simplest way to accomplish this is to modify the shift operator S, equation ([3|) 
while keeping the coin operation as before, equation ([2]). The variable shift operator S{x,t) is then 
composed of the following three variants in addition to the regular shift in equation ([3]) that is 
applied at positions with both edges present. If the edge in the plus direction is missing, 

S+| — 1, x) = \ — l,x — 1) 

S+|+l,x) = |-l,x), (8) 
if the edge in the minus direction is missing, 

S_ | — 1, x) = 1+1, x) 

S_|+l,x) = |+l,x + l), (9) 

and if both edges are missing, 

So|-l,x) = 1+1, x) 

So|+l,x) = |-l,x). (10) 
Since S+ a nd S_ occur in pa irs either side of gaps, the overall S(x, t) is unitary. This model was first 



studied by lRomanelli et al.l |26l ]. who found that, as the proportion of gaps increases {p decreases), 
the distribution of the quantum walk tends to the classical binomial spread, after a characteristic 
time that scales as 1/(1 —p). This is illustrated in figure [21 (which has been plotted for single runs, 
i.e., no averaging over the different possible sequences of missing gaps). Dynamic gaps thus act 
as a type of decoherence that transforms the quantum walk behaviour into classical random walk 
behaviour. This is thus an example of behaviour with a single scaling of ^/t, but a prefactor that 



varies depending on the parameters (the percolation probability in this case). iRomanelli et al.j [26] 
compare this with decoherence due to periodic and random measurements (non-unitary evolution) 
applied to both coin and position, w hich also reduce s the scaling to \/t but with larger prefactors 
Similar behaviour was also found bv iB run et al 



driving a quantum walk on the line. .Brun et al 



[33, l39( for decoherence applied only to the coin 
solved their model analytically for the spreading 
rate. They found that the limiting distribution scaled classically as but the prefactor was larger 
and determined by the decoherence rate. 



2.2 Quantum Tunneling 

Tunneling is a generic property of quantum systems. One of the easiest demonstrations uses light 
undergoing total internal reflection in a block of glass. If another block of glass is brought close 
to the surface where the reflection is taking place, some of the light can be seen to be transmitted 
even when there is a small air gap between the blocks. The probability of finding photons outside 
the glass at the point of reflection isn't zero, it falls away exponentially from the interface (the 
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Figure 2: Probability distribution after 10,000 time steps for a quantum walk on the line with 
dynamic gaps where the probability of an edge being present is p = 1 (black), 0.75 (red), 0.5 
(green), 0.25 (blue) and 0.0 (orange), using a Hadamard coin, equation ([2]), and a symmetric initial 
state (1+1,0) +i|-l,0))/^/2. Only even positions shown since odd positions are unoccupied. 



evanescent wave). If a second glass block is placed close to the first one, the photons return to 
propagating in a normal wave-like fashion in the second block. The amplitude of the transmitted 
wave is given by the amplitude of the evanescent wave at the second interface. 

We can make a simple model of quantum tunneling by allowing the quantum walker to hop 
over a broken link using the biased coin C2'''* from equation ([7]). Changing rj is equivalent to 
changing the size of the air gap between the glass blocks. The quantum walk dynamics will now 

have a position dependent coin operator, while the shift operator remains as in equation ([3]). We 

(v) 

apply C2 at the sites on both sides of the missing edge to create the tunnel. Using < 77 < 0.5 
corresponds to tunneling with a reduced probability of transmission, while 0.5 < r] < 1 would give 
us enhanced transmission, which we don't have a use for in this model. A coin operator composed 
of a combination of these two operators, c'^^'^^ and is easily seen to be unitary, because the 
full coin operator is a block diagonal matrix, C (8> U^i and the position dependence simply inserts 
the different C into the appropriate blocks. 

To see how the tunnelling works in more detail, suppose there is a gap between position x and 
position X -\- 1. If the walker arrives at x from x — 1, the coin will necessarily be in state |+1). 
Applying a step of the walk, \J{r]) = SCg''^ to the state |+l,x) gives 

S&^''\+l,x) = Vl-r/|-l,x-l) + ^|+l,x + l). (11) 

A fraction of the amplitude has thus crossed the gap. We now need to apply the same operator 
U(?7) to this amplitude at position x + 1, 

SC^'')^|+l,x + l) = Vr/(l-7?)|-l,x) +7?|+l,x + 2). (12) 

Thus after two steps, only a fraction rj of the amplitude has progressed completely past the gap. 
This compares with the expected fraction of half the amplitude that would have reached this 
position after two steps of the walk if the gap weren't there. Rigorously defining the probability of 
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tunnelling through the gap requires a more thorough treatment of the delocalized quantum walker, 
for example, using a traveling wave incident on the gap, rath er than localised position states. The 
coin operators then act like beam splitters, see lKnight et alj 3l|. For our purposes we don't need 
a precise definition, the biased coin gives us a simple way to interpolate smoothly between an 
impassable gap and the unobstructed walk, which allows us to explore the intermediate behaviour 
qualitatively. By placing a pair of gaps at ±20 and starting the quantum walk at the origin, 
the reflection and transmission at the gaps is nicely illustrated in figure [3l A different model of 




^ "^'tice units) 



-100 



100 



Figure 3: Probability distribution for a ID walker on a line showing the time evolution over 100 
steps using a Hadamard coin, equation ([2j), and a symmetric initial state (|-|-1,0) -|- i| — 1,0))/\/2. 
The tunneling strength rj = 0.25, and gaps flve edges wide are placed at ±20, as indicated by the 
labeled mark at t = 0. Only positions with non-zero probability of occupation are shown, since 
odd positions are unoccupied at even time steps and vice versa. 



tunneling to account fo r imperfect shift operations in a qua ntum walk on an uri broken line has been 
studied numerically by lDiir et al.l 40(], and analytically by lAnnabestani et al.l 411 ] . who concluded 
that in this scenario, the extra shifting does not destroy the quantum spreading rate. 

Figure [H shows the effect of tunneling on a quantum walk on a line with randomly placed static 
gaps and a tunneling strength of r/ = 0.25 (half the normal unbiased coin transmission rate). The 
number of missing edges is controlled by the probability p of any given edge being present. After 
100 steps, the spreading is essentially still at the fully quantum rate for all values of p. By 1000 steps 
a small central peak has appeared at the expense of the leading edges, while for 10,000 steps the 
intermediate values of p have reduced the walk to classical-like distributions (like classical random 
walks on the unbroken line). The random distribution of gaps in the percolation lattice removes a 
random subset of the possible same-length paths between any two points. This slowly destroys the 
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Position on line 

Figure 4: Probability distribution for a ID walker on a line with static gaps and a tunneling 
strength of = 0.25 after 100 steps (top), 1,000 steps (middle), and 10,000 steps (bottom) using 
a Hadamard coin, equation ([2]), and a symmetric initial state (|+1,0) + z| — 1, 0))/\/2, where the 
probability of an edge being present is p = 1 (black), 0.75 (red), 0.5 (green), 0.25 (blue) and 0.0 
(orange). Only even positions shown since odd positions are unoccupied. 



• p= 1.0 
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Figure 5: Probability distribution for a ID walker on a line with dynamic gaps where the probability 
of an edge being present is p = 1 (black), 0.75 (red), 0.5 (green), 0.25 (blue) and 0.0 (orange), and 
a tunneling strength of r] = 0.25 after 10,000 steps using a Hadamard coin, equation and a 
symmetric initial state (|+1,0) + i\ — l,0))/^/2. Only even positions shown since odd positions are 
unoccupied. 



long-term coherence required for the full quantum spreading rate. After 10,000 steps, only p = 1 
(no gaps) and p = (all tunneling) show double-peaked quantum speed up. Figure [5] shows the 
combination of the same dynamic gaps as in figure [2] and tunneling with rj = 0.25. The tunneling 
clearly enhances the progress of the walk considerably, but the dynamic gaps still cause decoherence 
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and a consequent return to a classical (binomial) shaped distribution for the intermediate values 
of p. However, the width of the distribution is significantly larger than the classical random walk 
distribution on the full line. The standard deviation for the classical random walk for 10,000 steps 
is 100, while the standard deviation for the quantum walk with tunnelling strength p = 0.25 is 479, 
for p = 0.5 it is 397, and for p = 0.75 it is 414. Tunnelling is thus allowing the quantum coherences 
to survive for longer. 

In contrast to the decohering effects of random gaps, using a regular pattern of static gaps at 
every other edge (like a diffraction grating) slows the progress of the walk, but does not affect its 
characteristic quantum shape, see figure [6l This is an illustration of the sensitivity of a quantum 




Figure 6: Probability distribution for a ID walker on a line showing evolution over 100 steps using a 
Hadamard coin, equation ([2]), and a symmetric initial state (1+1,0) + i| — 1, 0))/\/2. The tunneling 
strength rj = 0.25. Gaps are placed at every other link. 



walk to symmetry. Only half that density of random gaps significantly disrupt s the quantum walk, 
as shown in figure HI Regular patterns of varying coin bias have been studied by lLinden and Sharam 
(43 |. who find quantum spreading in most cases, b ut with widely varying spreading rates. General- 
isation of this model by lShikano and Katsural 431 ] suggests that localization vs spreading depends 
on whether the variation in the coin is governed by a rational or irrational number. 

In order to make quantitative comparisons of the effects of varying the percolation parameter 
p and the tunneling strength we calculated the root mean squared (rms) value of x for the 
probability distribution, i.e., Xj-ms = Y^a;2). We will use the overbar to denote averaging over a 
single quantum walk. This gives us a measure of how fast the quantum walk is spreading from its 
starting point at the origin x = Q. We could equally well have used but Xrms (often referred to 
as the standard deviation when the mean is zero) is more commonly used for this purpose. Since 
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we want to know about the average behaviour, we repeated the quantum walk for 1000 different 
random distributions of gaps, but for a walk of only 40 steps, to ensure we obtained reasonable 
statistics. We thus have (a^rms) = (v^^^)) '^^^ measure of the spreading of a quantum walk on a 
line with random gaps. We will use (.) to denote averaging over many different randomly generated 
percolation samples. It is also useful to examine how variable the spreading is for fixed p and rj, 
so we also calculated (ri^XYms 

), the standard deviation of Xrms over the set of random percolation 

samples. 

For a random arrangement of missing edges that stays constant throughout the walk, (xrms) and 
(^{xrms) are shown in figures [7)^ a) and[71^b). The (xrms) in figure [7|^ a) is as one would expect, with 
the gaps having a big impact on the spreading for low tunneling rates, while even a small amount 
of tunneling allows the spreading rate to steadily recover. The cT(xrms) in figure W[h) shows that 
this behaviour is reliable (low variation) except where both the number of gaps and the tunneling 
rate are low. 

For dynamic gaps, the equivalent pair of plots are shown in figures [7|^c) and[7|^d). The in- 
troduction of dynamic gaps can be seen to have a steep but less drastic effect than static gaps 
by comparing figure [TJa) with figure [7)^c). The (T(xrms) is also lower, indicating the behaviour is 
less variable between different lattices with the same percolation probability p. This is due to the 
dynamic gaps providing an averaging effect during the course of the walk. 

The main message from these results for quantum walks on a line with edges missing is that a 
random pattern of missing edges leads to decoherence in the quantum walk, reducing the spreading 
to the classical random walk scaling proportional to \/t, but with varying prefactors. However, the 
crossover to classical scaling hap pens quite slowly, after a crossover time proportional to 1/(1 —p) 



as found bv lRomanelli et al.l [32l | for the case of dynamic gaps. This means the decoherence effects 
are minimal for smaller systems allowing the full quantum speed up to be exploited, while larger 
systems will only benefit from a prefactor enhancing the classical scaling. 



3 Walk on a two-dimensional grid 

A two-dimensional lattice allows a wider range of behaviour to be studied, without resorting to 
quantum tunneling. Since there is more than one path between two points on the lattice, the walker 
can take a (probably longer) path around the gaps when they are at a low density. Each vertex in 
a 2D lattice has four edges connected to it, so the quantum walker now needs a four-dimensional 
coin. We label the four directions \L), \R) corresponding to | — 1) and | -|- 1) on the x axis and \D), 
\U) for I — 1) and | + 1) on the y axis. The operator for "tossing" the coin can now be any 4x4 
unitary. The most common choice is based on Grover's diffusion operator. 
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The general form of Grover's diffusion operator for a d-edged vertex is 2/d — U^, it happens to 
be unbiased only for d = 4. It is chosen because of its symmetry, the incoming edge receives a 
sign fiip, but the other three edges are treated equally. There are many possible unbiased coin 
operators based on generalised Hadamard matrices 4J] , but these will have other phases that treat 
the outgoing directions differently. 
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(c) {xrms) for dynamic gaps. . r . 

(d) a{Xrms) lor dynamic gaps. 

Figure 7: Plots of (xmis) (left) and o"(xrms) (right) against < t] < 0.5 and < p < 1.0 for a ID 
quantum walk, with static (top) and dynamic (bottom), run for 40 steps, and averaged over 1000 
random gap arrangements. Note scale mismatches in vertical axes. 

To complement this choice of coin operator, the shift operator S4 must also invert the coin 
states, 

Si\L,x,y) = \R,x- l,y) 
S4\R,x,y) = \L,x + l,y) 
Si\D,x,y) = \U,x,y - 1) 

S4\U,x,y) = \D,x,y + l). (14) 



12 



This ensures that the coin operator is then apphed to a coin state that labels the edge from which 
the walker arrived, which is the one that receives the minus sign from the Grover coin operator. 
For example, 



cf''"^\L,x,y) = ^ {-\L,x,y) + \R,x,y) + \D,x,y) + \U,x,y)) . (15) 

First we recall the behaviour of a quantum walk on a lattice with all edges present. This is 
illustrated in figure [8l For most choices of coin and initial state, the quantum walk tends to stay 




Figure 8: Probability distribution for quantum walk on 2D lattice run for 40 steps with initial 
states 0.5(|L,0) + |i?,0) + |C/,0) + ID, 0) (left) and 0.5(-|L, 0) - 0) + |C/, 0) + |Z), 0) (right). Only 
even positions are shown since odd positions are unoccupied. Note scale mismatch in vertical axes. 



around the origin. For a well-behaved coin (with a good degree of symmetry), there is one precisely 
chosen state for which it spreads out in a ring of variable height that resembles a circular version 
of the walk on the line [Sj]. If we start the walker off from the origin (x = 0, y = 0), we can use the 
average distance of the walker from the origin f, the average of r = a/x^ + y^, as a measure of the 
spreading. This is playing the same role as Xrms for the walk on the line, although f is more directly 
comparable with the average of since the square root is taken before the averaging. In figure 
IHl the cases with the two extreme values of f for the Grover coin operator are shown, on the left f 
after 40 steps is 33.1, while on the right it is only 7.2. Nonetheless, this still represents a quantum 
spreading rate, despite the apparent dominance of the central peak. The classical spreading of ^/t 
is only 6.3 for 40 steps. Other initial states show a mixture of the ring and central peak, with 
corresponding intermediate spreading rates. However, the central peak is present and dominates 
over the ring for all except the symmetric starting state. This is in sharp contrast with the walk 
on the line where the spreading rate is the same for all initial coin states, only the symmetry of 
the distribu tion changes. A i ialy tic t reatment of the qu antum walk in two dimensions has been 



provided by iGrimmett et al.l [45| and iGottlieb et al, 
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3.1 Two-dimensional percolation lattices 

Applying a quantum walk to a two-dimensional percolation lattice is a straightforward extension of 
the techniques already described for varying the dynamics of the walk on the line. For both bond 
and site percolation, the resulting lattice will in general contains vertices of degree 4, 3, 2, 1 and 
0. Each needs an appropriate coin operator. We constructed these in the most symmetric possible 
way, basing the lower dimensional coins on the lower dimensional Grover operators padded with 
(minus) the identity for the directions where the edges are missing. For example, if the \U) edge is 
missing, the coin operator becomes 



C. 



/-I 
2 

2 



while if both the \U) and \L) edges are missing, 



2 




2 \ 

2 

-1 

-3/ 



(16) 



(LU) 



















1 








1 








vo 








-V 



(17) 



Provided the initial state is correctly chosen, so the subspace of the missing edges is never populated, 
it doesn't matter what appears in those rows and columns of the coin operators. However, if one 
wishes to interpolate between missing edges and the full lattice, in a generalisation of the tunnelling 
studied on the line in the previous section, the sign of the entries on the diagonal must be the same 
throughout the coin operator, as shown above. For consistency, the coin operator for all but one 
edge missing 

/-I \ 
^{id) ^0-10 
4 0-10 

V -1/ 

applies a minus sign to the amplitude before it returns along the edge. 

Since the initial coin state has a significant effect on the spreading rate of the walk, we therefore 
tested several different initial coin states, to see how this affected the behaviour, in combination 
with the missing edges in the percolation lattices. We tested the initial states iV'max) and iV'min) 
that give the maximum and minimum spreading rates as shown in figure [HI 



(18) 



0.5(|L,0) + |ii,0) + |C/,0) + |i?,0), 



(19) 



and 



0.5(-|L,0) - |i?,0) + |C/,0) + |Z),0). 



(20) 



These differ only in the phases between the different coin directions. We thus also tested random 
phases of ±1 between the different coin directions. 



0.5(±|L,0) ± |i?,0) ± |C/,0) ± |D,0). 



(21) 
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This is a similar strategy to that used by iTregenna et alj 341 ] when obtaining the maximum and 



minimum spreading rates. Restricting the variation to phases avoids skewing the walk, which would 
systematically increase the average distance traveled, making it harder to extract the generic average 
behaviour from the results. In all cases, the walker started at the origin. If there happened to be 
missing edges, those coin states were deleted from the initial state, and the remainder renormalized. 
If all the edges happened to be missing around the origin, a spreading of f = was recorded, since 
no walk was possible. 

For each initial state in equations (]19p to (j2ip . we performed a numerical simulation of a 
quantum walk on both site and bond percolation lattices, averaging over 5000 randomly generated 
lattices of each type. We were able to run up to 140 time steps on a 3GHz processor in reasonable 
time (a week per 5000 sample lattices per initial state). With (2 x 140 + 1)^ = 78961 lattice sites, 
this is equivalent in memory requirements to nearly 40000 steps of a walk on a line. Arguably, 
5000 sample lattices is not a sufficiently large sample for such large lattices. We investigated the 
sensitivity to sample size by comparing averages for smaller lattices (40 to 80 time steps) for 1000, 
5000 and 20000 samples. The larger samples produced smoother averages as expected, but the 
smaller samples gave good enough results to draw conclusions from. 

Figures [H [10] and [TT] present the results of these simulations in three complementary ways. 
In figure [9l we show the results for quantum walks with random phases in their initial states. 
The average distance from the starting position at the origin for a single quantum walk, f where 
r = Y^a;2 + y^, is given by 

f = Y,.Jx^ + y^\i^{x,y)\\ (22) 
x,y 

Here \il){x,y)\'^ is the probability of finding the particle at the site (x,y), i.e., summed over all 
coin states, and the sum over x,y is over all lattice sites that could be reached by the walker, i.e., 
—t<x,y<t. Clearly, f can vary, depending on the exact arrangement of missing sites or edges in 
the percolation lattice. We have therefore further averaged over 5000 random instances of bond or 
site percolation lattices, and plotted the resulting (f) as a function of the percolation probability 
p, for different lengths of quantum walk from 10 to 140 steps. Values of (f) below one that occur 
for small p are not interesting: the low density of sites or edges in the lattice means that most of 
the time the quantum walk is unable to move from its starting state. While there is clearly some 
noise in the data, particularly for site percolation, the trends are clear. At around p = 0.4 for bond 
{p = 0.5 for site) percolation, (f) starts to rise. This is still below the critical probability pc = 0.5 
{pc = 0.59 ... ) and the value of (f) is not sensitive to the number of time steps because in general 
it is still confined to a small local region and cannot see the full size of the lattice. We then observe 
the lines for different time steps separating out below the main bundle in succession, smallest first, 
and flattening off towards the corresponding value for the complete lattice at p = 1. This is a nice 
illustration of finite size effects. The size of the localized clusters grows as p approaches pc from 
below. As the cluster size becomes comparable with the size of the walk (10, 20 time steps, etc.), 
the limitation of (f) due to cluster size is replaced by limitation due to the number of steps, and 
the walk is no longer sensitive to increasing cluster size. Walks of about 100 steps or longer are still 
close to synchronized at p = pc, indicating that quantum walks of this size should provide good 
estimates of critical parameters. 

In the region where pc < p < 1, the rate of growth of (f) slows then steepens again. For 
Pc ^ 0.95, the walk no longer distinguishes between bond and site percolation, (f) is the same 
for both. This is somewhat surprising, because removing a site from the lattice removes four 
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Figure 9: Comparison of the average distance (f) from the starting point for bond (upper set) and 
site (lower set) percolation for the number of steps from 10 to 140 as shown in the key, as a function 
of the percolation probability p, with initial states using random phases, equation ([2T|) . 

neighbouring edges, compared to removing a single edge for bond percolation. It suggests that 
in the regime where missing sites or edges can be expected to be isolated from each other, they 
both behave as lattice defects, with the broken lattice symmetry more important that the detailed 
structure of the defect. 

To understand this behaviour in more detail, we turn to figure [TOl where walks of 60 time 
steps are compared for all three initial states. Also plotted is cr(f), the standard deviation of f, to 
quantify how variable the spreading is between different percolation lattices. Note first that the 
average spreading (f) for p = 1 ranges from around 14 for iV'min) to 35 for IV'max)- The value for 
IV'ran) IS around 18, half that for IV'max)- However, a{f) for iV'ran) is around seven, indicating that 
the average value conceals significant variation over the random phases. This is as expected, since 
the randomly chosen phases should include both iV'min) and iV'max)- The position of \^pra.n) close to 
IV'min) indicates the weight of the distribution lies nearer the minimum spreading rather than the 
maximum. All this is already known from previous studies, of course, but it is useful to recall the 
details before considering the behaviour for p < 1. 

Looking at the standard deviations, it can be seen that the variation for iV'ran) is larger for site 
percolation, typically around five compared with around two for bond percolation for pc < p < 0.95. 
Thus, although the average behaviour is similar, missing sites produce more extreme cases (high 
and low) for the same number of missing bonds in bond percolation. This suggests further work to 
elucidate the source of this difference would be fruitful. 

The variation of (f) over time with p can be plotted in 3D, see figure [TT] for two examples. This 
shows how similar the behaviour is for both bond and site percolation, with any choice of initial 
state, but it does not bring out the most interesting features to do with how the scaling varies with 
the percolation probability. The necessary scaling analysis is described in the following subsection. 
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Figure 10: Comparison of the mean distance (f) from the starting point for bond (left) and site 
(right) percolation over 60 steps, averages over 5000 random percolation lattices, as a function of 
the percolation probability p, with initial states giving the maximum and minimum distance in the 
full lattice case, plus averaging over random phases in the initial state. The standard deviation 
(T(f) over the random percolation lattices is also shown. 




Figure 11: Variation of spreading rate (f) with percolation probability up to 140 steps for bond 
(left) with max spread initial state, and site (right) percolation on 2D Cartesian lattice with random 
initial state. Note difference in vertical scales. 

3.2 Scaling analysis 

A log-log plot of the same data in figure [11] (right) is shown in 2D in figure [121 Apart from a 
regular zigzag over a cycle of four steps, due to the lattice and coin symmetries, the line for a fixed 
value of p is very close to a straight line, with slope varying from zero (small p) to one (p = 1). 
This suggests the data will be well fit by (f) oc with < a < 1. Mindful of the finite size effects 
discussed in relation to figure [H fits to determine a were done using only the data for time steps 
100 to 140. The resulting values for the exponent a are shown in figure [TSj The values of a agree 
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Figure 12: Variation of spreading rate (f) with percolation probability up to 140 steps for site 
percolation on 2D Cartesian lattice with random initial state. A log-log plot so slope gives a 
where (f) ~ t". 



fairly well with the expectation that the first rise above zero should be seen around the critical 
point pc = 0.5 (bond) and pc = 0.59 . . . (site). The initial state that produces maximum spreading 
has a noticeably higher exponent, around 0.05 larger than for random phases or the minimum 
spreading state, from the critical point through to near p = 1, implying that the advantage of the 
optimal initial state is not just a prefactor on the same scaling rate. All initial states follow the 
same overall pattern of a steady rise from the critical point, a flattening off towards a = 0.5, the 
classical exponent for the perfect lattice, in the region of p = 0.9, followed by a steep rise to the 
quantum spreading rate of a = 1. 

To test the influence of finite size effects on these results, the same fits were done using more 
of the data, with t = 80 and t = 60 as lower cut off points in place of t = 100. The slightly higher 
exponents these fits produce are shown as an inset in figure [T3l along with more detailed data, for 
the case with random phase initial states. These show that the trend towards lower values of a 
is slow, and hence, based on these simulations, we cannot make a strong prediction for what the 
asymptotic variation of the exponent with percolation rate is going to be. Nonetheless, we can 
confidently say that for systems of limited size, there is a significant window for p > 0.9 in which 
the quantum advantage exists in the form of faster spreading than a classical random walk. 



4 Summary and discussion 

Even for the simple case of a coined quantum walk on the line we find contrasting effects brought 
about by the interplay between missing edges and quantum tunnelling, in the presence of ordered 
or disordered arrangements of the missing edges. Ultimately, for long times, the disorder of random 
missing edges causes the spreading to reduce to the classical ^/t rate. Yet for shorter walks, even up 
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Figure 13: Fractional scaling exponent for 2D percolation lattices derived from from data for 
t = 100 to 140 steps. Initial states (bond, site) for maximum (blue, cyan), minimum (red, 
orange) spreading and random phases (green, lime). Inset shows greater detail for 0.8 < p < 1.0 
for random phases, and results of using lower cut off of 80 (dashed) or 60 (dotted) steps to 
determine exponent. 



to a thousand steps, the spreading is still largely dominated by quantum effects. And for ordered 
arrangements of missing edges, the quantum spreading rate is maintained, though with varying 
prefactors. These res ults are consisten t with those of a number of other analytical and numerical 



stu dies, in particular lRoma,nelli et al.l 32l | who first studied randomly changing gaps on the line. 



and lLinden and SharamI |42l ]. who provide detailed analytic treatment of periodic ally varying coins . 



The connections between disorder and decoherence were discussed in detail by iRomanelli et al. 



highlighting the role of randomness from whatever source in degrading the quantum coherence. 
The scaling of quantum walks on the line with decoherence or other sources of disorder are in 
general either linear spreading or \/t spreading, with different prefactors depending on the details 
of the dynamics. Intermediate scaling in one-dimensional walks is not obtained except transiently 
as linear crosses over to a long time limit of y/i. 

In contrast to the ID cases, on 2D percolation lattices we find the quantum walks show frac- 
tional scaling of the spreading, i.e., f oc for pc < p < 1, with < a < 1. Classical random 
walks on percolation lattices also exhibit fractional scaling behaviour, with < a < 0.5. One 
interesting question is whether this behaviour seen for finite-sized walks extends to the long time 
limit. Although our simulations seem to have reached a size where finite size effects are small, the 
slow convergence, coupled with the lessons from one dimensional quantum walks where decoherence 
dominated only after thousands of steps, means we do not think we can confidently predict the 
large t behaviour. For example, a possibility consistent with our results is that the steep rise in 
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the region of 0.95 < p < 1.0 will becomes a "step" function at p = 1 as t — )• oo, and that the 
portion below a = 0.5 will follow the scaling for classical random walks. The intuition for this is 
that the randomness in the percolation lattice would thus again act as decoherence in the large 
t limit. Whether this is correct depends on whether a very small density of defects is enough to 
disrupt the entire walk in the whole two-dimensional lattice in the long time limit. Studies of 
absorbing boundaries and trapping are relevant to this question, in one dimension the presence of 
traps does not inevitably catch a quantum walker [33], while a classical walker will eventually fall 
into a single trap in both one and two dimensions. This will be an interesting question to answer be- 
cause trapping behaviour is particularly important for models of exciton transport, as discussed by 
47|. an d trappi ng can also be exploited to guide the exciton to the desired location. 



Miilken et al 



Feldman and Hilleryl 481 ] and iHillery et al.l 491 ] provide a scattering matrix method for treating 
finite graphs analytically, and determining solutions that traverse the graph or are trapped within 
it. This may provide a route to understanding how structures in the percolation lattice contribute 
to trapping the quantum walker. 




ran phase 
ran phase 

classical 
classical 
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Figure 14: Fractional scaling exponent for 2D percolation lattices derived from from data for 
t = 100 to 140 steps. Initial states (bond, site) for random phases (green, lime) are shown (same 
data as in figure [T3l) . Classical random walk analysed in the same manner is shown for comparison 
(maroon, magenta), which displays obvious finite size effects. 

Aside from the question of trapping behaviour, we offer some further intriguing observations 
that suggest directions for future study. Figure [U] shows the data for random initial states from 
figure [13] plotted along with numerical results for running a classical random walk on percolation 
lattices up to 140 steps, and processed to extract the scaling in exactly the same way. The first 
and most obvious problem is that 140 steps of a classical random walk still suffers from large finite 
size effects! This is easily understood, f for a classical random walk on the fully connected lattice 
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is ^/t, which is about 12 for t = 140. This is more comparable with t = 40 for a quantum walk 
(with random phases in the initial state), as can be obtained from figure El reading off the values 
for p = 1. Since f gives an indication of the size of the largest structures that the quantum or 
random walk will be sensitive to, if the missing sites or edges are on average further apart than 
that, the walker will hardly notice any before finishing the walk and reporting its position. Working 
in the other direction, a quantum walk of 100 steps (with random phase initial state) has f = 30. 
If we want a classical random walk with f = 30, will need to run it for 30 x 30 = 900 steps! The 
advantage of the linear scaling with t of the quantum walks is suddenly very apparent. 

As already noted, we appear to have achieved a reasonable convergence of finite size effects 
on percolation lattices by running a classical simulation of a quantum walk for 100 or more steps. 
These simulations required around a week of processor time on a fast workstation, a fairly modest 
amount of computational resources. This suggests that using such simulations to obtain numerical 
values for percolation lattice parameters is worth investigating. In other words, this could provide 
a route to improved classical simul ation methods. There are other examples of quantum algorithms 
providing better classical methods. iLove and BoghosianI 501 ] discuss a number of examples in detail. 

However, the main reason for running the classical random walk was to compare with the 
quantum walk and further understand our results. Much larger classical simulations were not 
possible with our available computational resources, we would have needed an order of magnitude 
more computational power, or to implement more efficient simulation techniques. Comparisons 
using the results in figure O can only be suggestive, nonetheless, the suggestions are intriguing. 
If we assume that the shape of the classical curves in figure [14] won't change much, so they will 
simply shift to the right as finite size effects reduce, then we can compare the shape with the 
curves for quantum walks. For site percolation, the classical and quantum curves are the same 
shape up to about p = 0.9 on the quantum curve. The corresponding curves for bond percolation 
are not, however, the quantum curve lies significantly below the classical curve for the region 
0.55 < p < 0.85. This adds further evidence for basic differences in quantum walk behaviour 
on bond and site percolation lattices that the differences in standard deviation noted in figure [TOl 
imply. In particular, that quantum walks are slower than classical in some regions, whi le faster than 



classical in others. Such behaviour has been noted in other contexts, see for example I Agliari et al 



51|. 

Nonetheless, the main message we want to emphasize from this study is that for modest system 
sizes of a few hundred sites, and defect densities below 10%, faster-than-classical fractional scaling is 
very much the dominant feature. This is an interesting regime with many applications in both bio- 
and nano-scale materials, where quantum effects are important to fully understand their properties 
and behaviour. 
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